Minimal conditions for local pure-state entanglement manipulation 
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We find a minimal set of necessary and sufficient conditions for the existence of a local procedure 
that converts a finite pure state into one of a set of possible final states. This result provides a 
powerful method for obtaining optimal local entanglement manipulation protocols for pure initial 
states. As an example, we determine analytically the optimal distillable entanglement for arbitrary 
finite pure states. We also construct an explicit protocol achieving this bound. 
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The existence of non-local correlations, or entangle- 
ment, between parts of a composite quantum system is 
at the heart of quantum information theory and its ap- 
plications Q. In recent years, much effort has been ex- 
pended on the problem of how to define and quantify the 
entanglement of a given state in physically meaningful 
ways. One very fruitful approach, first pursued by Ben- 
nett and co-workers Q-^], is to regard entanglement in 
terms of the limitations that exist to the manipulation of 
a composite system when each subsystem is operated on 
locally. A paradigmatic situation is as follows: suppose 
Alice and Bob each possess part of a quantum system, 
which is prepared in a state p. Qualitatively, the exis- 
tence of entanglement implies that some transformations 
of p which arc in principle possible cannot be realised 
if Alice and Bob are only allowed to perform local oper- 
ations on their respective subsystems, and to exchange 
classical communication. (Transformations of this type 
can be referred to as 'local transformations', or 'LQCC 
for short). 

A quantitative way of expressing this fact is in terms of 
so-called entanglement monotones (EMs) ||. These are 
functions e (p) of the quantum state that can, on average, 
never increase under LQCC There are many known 
EMs, for example the entanglements of distillation (^-^] 
and formation , and the relative entropy of entangle- 
ment Q (in fact, any reasonable measure of entanglement 
must by definition be an entanglement monotone, and 
vice- versa). Despite their different physical interpreta- 
tions, they all share a common feature: a transformation 
which, on average, increases any single EM cannot be 
realised locally. In other words, they provide necessary 
conditions any local transformation T must satisfy. 

A natural question that presents itself is then: what 
are sufficient conditions for T to be local? In other words, 
we would like to have a set {e^} of entanglement mono- 
tones such that, if the average (ei(T[p])) < £i(p) for 
all i, then T is local. Ideally, this set should also be 
minimal, in the sense that these conditions should not 
be redundant JjJ. An important result in this direction 
was recently presented by Nielsen | px| , who found suffi- 
cient conditions for the locality of transformations that 



take one given pure state to another with 100% probabil- 
ity. In the present Letter, we extend Nielsen's theorem 
to the case where the transformation need not be deter- 
ministic, that is, when T may lead to several possible 
final states. We demonstrate that, for this case, a set 
of EMs recently introduced by Vidal || is in fact mini- 
mal in the sense described above. They therefore provide 
us with a powerful universal tool for finding optimal lo- 
cal entanglement manipulation protocols. We apply it 
to the problem of entanglement concentration (or purifi- 
cation), which concerns understanding to which extent 
distant parties can extract a maximally entangled state 
from a non-maximally entangled one using only LQCC 
[p|-fI],[rT| . This is one of the central problems of quantum 
information theory, and is of crucial importance for all 



applications, such as teleportation |12|, which require the 
existence of maximally entangled states between distant 
parties. With the help of our generalisation of Nielsen's 
theorem, and of results from the well-known simplex op- 
timization method of linear programming theory fl3f| , we 
are able to analytically determine the optimal purification 
protocol for the case where Alice and Bob share a given 
pure state By 'optimal' we mean the following: as- 
sume that Alice and Bob locally manipulate their shared 
state until they obtain cither a maximally entangled state 
(of some dimension), or a completely disentangled one. 
We determine the strategy that awards them, on average, 
with the largest amount of distilled entanglement, which 
we find to be 
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where ol\ > ... > ojjv are the nonzero Schmidt coefficients 
of|# ' 

It is important to stress that our results pertain to 
any finite shared state. Until now (see also endnote), the 
problem of finding the best purification protocol in the 
sense above had been completely solved (for pure states 
and some particular mixed states), only in the asymp- 
totic limit, where Alice and Bob share N — > oo identical 
copies of the same state |^,|| . This limit has fundamental 
significance in quantum information and communication 
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theory, for example for deriving bounds on channel ca- 
pacities 0. Nevertheless, it is clear that in practise Alice 
and Bob will always share only a finite, in general small, 
amount of entanglement. Thus, as a number of authors 
P,p| JlO|JI^ | have stressed, it is also important to under- 
stand entanglement transformations in this regime, with 
the asymptotic results emerging in the suitable limit. 

Suppose then that Alice and Bob share a pure state \ip) 
of a bipartite quantum system, with ordered Schmidt de- 
composition \ip) — YliLi V®i\* A ) ll^l- Vidal |5] has 
shown that each of the following set of functions of the 
on constitutes an entanglement monotone: 



N 



Ei(\i ) )) = ^2a i , 1<1<N. 



(2) 
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We can use these monotones to describe the follow- 
ing theorem due to Nielsen fiofl : let |?7) be another 
pure bipartite state. Then there exists a local trans- 
formation that takes \tp) to |?7) with 100% certainty iff 
Ei(\i])) < Ei(\ip}), 2 < I < N. In other words, the 
{Ei} form a sufficient set of monotones for this kind of 
transformation. In fact, since they also uniquely deter- 
mine the Schmidt components of \xf>) and |?7), which com- 
pletely and minimally characterise such transformations 
(sec. 5.1 of ||), it follows that {E{\ is actually a minimal 
set of EMs in this case. 

Quantum mechanics is not, however, concerned only 
with deterministic transformations. As long as Alice and 
Bob do not lose or discard information about their sys- 
tem, the most general transformation they can apply on 
\tp) will produce one of m possible pure states \r)i), with 
probability pi. We demonstrate now that Vidal's mono- 
tones also provide necessary and sufficient conditions for 
these general transformations to be realised locally. 

Theorem 1: Let 2 distant parties share a pure state 
\ip) ; let {Ivj}}™^ be a set of m other pure bipartite 
states. Then a transformation T of that outputs state 
\r/j) with probability pj {^jPj = lj can be realised us- 
ing LQCC iff the N entanglement monotones Ei do not 
increase on average, that is, iff 



J2Pj E d\Vj))<Ei(m,l<l<N. 
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Proof: Necessity follows from the definition of an en- 
tanglement monotone, and is proven for the Ei functions 
in i{. To prove sufficiency, assume eq. (||) is satisfied. 
We will construct an explicit local strategy that realises 
the transformation T. First of all, it is clear that we only 
need to consider the special case where all target states 
\rjj) have the same Schmidt basis as (which we can 
refer to as the 'standard' basis). The general case then 
follows from the following simple facts: (i) any two states 
with the same Schmidt components are interconvertible 



by a local unitary operation, so that to realise T one only 
needs to generate, with probability pj, a state \r/j) with 
the same Schmidt coefficients as \rjj) in the standard ba- 
sis, (ii) If two or more target states , \r)j n ) have 
exactly the same Schmidt components, one can generate 
the state with probability Y^k=iPjki anc ^ then 'roll 
a classical die' with relative probabilities ■<- ^^ H to de- 

cide which one of the \f]j k ) to transform to. 

Suppose then that the target states can all be written 
in the ordered Schmidt form \r)j) = YliLi y/Pji Ka) |*b) 
(note that the number of nonzero Schmidt components 
of \rjj) cannot be greater than N JTJ]). 

Let us now define the average target state \fj) as 
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It can be seen that 7$ > 7i+i, so 

N m m 
i=l 3 = 1 3 = 1 

where we have used condition (||). We can therefore ap- 
ply Nielsen's theorem, which implies that there exists a 
local protocol L for deterministically converting from \tp) 
to 1 77). Let us now define the following set of positive 
operators on Alice's subspace: 
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\ia) (i A \ , 1 < j < m 
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We can see that, for 1 < j < m : 
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where P is the projector J2iLi \^a) Together with 

the complement 1a— P, the set {Aj}™_ 1 constitutes there- 
fore a local POVM which, if applied to \fj), outputs state 
\rjj) with probability pj. The combination of this POVM 
with the deterministic protocol L realises the required 
transformation T a 

This result can be directly extended to the case where 
the target states may be mixed. In this case, eq. (^) still 
holds (substituting pj for \rjj)), as long as we extend the 
definition of E\ using the 'convex roof rule || 

m 

Ei(pj)= min ^QijEiQrjij)) (9) 

ftf=2 J gyl'7<3>(»7«l j=i 



where the minimum is taken over all realisations of 
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Theorem 1 provides a powerful tool for optimizing lo- 
cal quantum transformations according to a wide range 
of criteria. For example, in ref. ||, the author seeks to 
determine the local transformation that maximizes the 
probability of converting one given pure state to another. 
He obtains an upper bound on this probability from the 
existence of the monotones Ei, and then constructs an 
explicit protocol realising the bound. Theorem 1 justi- 
fies this result, showing that a similar strategy will work 
for any optimization problem involving an initial pure 
state; in other words, the optimum given the constraints 
expressed in eq. (||) will always be achievable. 

We can immediately apply this result to the problem 
of optimally concentrating the entanglement of a finite 
bipartite pure state. This situation has already been con- 
sidered by Lo and Popescu [[u}, who have obtained the 
local protocol that gives the greatest probability of con- 
verting a given pure state to a maximally entangled 
state of any given number of levels. However, it may well 
be that Alice and Bob merely wish to concentrate their 
entanglement, without regard to what maximally entan- 
gled state they end up with. In this case, a reasonable 
question to ask is: out of all such local concentration pro- 
tocols, which one leads, on average, to the largest amount 
of shared distilled entanglement? 

The problem may be formally posed as follows: 
let Alice and Bob share a single pure state = 
\^ A ) \^b), whose entanglement they wish to 
concentrate using LQCC. Following the notation of ref. 
jLjJ l; let us define \<f>j) — Ei=i Ka) Kb) as a max- 
imally entangled state of j levels (note that \<f>i) is a 
product state). Consider the set of local transformations 
that generate \<j>j) , 1 < j < N, with probability pj [ p^[ . 
If we choose to measure the amount of entanglement in 
\4>j) by the von Neumann entropy oitrs \4>j) (<t>j I? namely 
Inj, then the average amount of distilled entanglement 
obtained from such a procedure is 

N 

(E) = Y, Pj ]nj. (10) 

7=1 

Our problem is to maximize this quantity over all prob- 
ability distributions for the pi that are consistent with 
the constraints in eq. (|3|). Theorem 1 then guarantees 
the existence of a local protocol leading to this optimal 
distribution. 

It is easily seen that, for I < j, 

E l (\4> j )) = J -^j±±, (11) 

and that it vanishes otherwise. In this case, the con- 
straints in eq. (§) read 

Ef=^(^)<Ef=^, i • I A • (12) 



This is a linear optimization problem with linear inequal- 
ity constraints, a kind widely studied in many fields of 
science and engineering. It can be solved using the tech- 
niques of linear programming theory, a branch of applied 
linear algebra that is familiar to most engineers, though 
not so well-known among physicists. We will not attempt 
to explain the terminology and results from this theory 
that are required for our solution; instead, we refer the 
reader to textbooks (e.g., Q). Our main result is 

Theorem 2: The optimal entanglement concentra- 
tion procedure for a single pure bipartite state \ip) with 
Schmidt coefficients a,\ > ...a/v > is one that pro- 
duces a maximally entangled state \<fij) of j < N levels 
with probability p° pt — j (<x, — aj +1 ). The correspond- 
ing optimal average distilled entanglement is (E) max = 

J2f=i ( a j ~ ot j+ i)jhxj. 

Proof: First, it is easy to check that this probability 
distribution satisfies (actually, saturates) all the inequali- 
ties in eq. ( |T2j ) . In matrix form, we have Bp = q, where p, 
q are vectors with components pj = j (aj — aj+i) ; qi = 
E f=i a j ' an d B is an upper triangular N x N matrix 
with components bij = for j >l, and otherwise. 

In the parlance of linear programming theory, this is a 
basic, feasible solution to the problem, with all the slack 
variables assuming the value zero. We can then apply 
the simplex algorithm to check whether this is the opti- 
mal solution or, if not, to find a better one. A sufficient 
condition for optimality ( jl3|], eqs. (2.36,2.37)) is that 
the following inequalities are all satisfied 

N 

z k = J2 CiPik > 0, l<fc<A^ (13) 
t=i 

where c$ = lnz is the coefficient of pi in eq. (|l0|), and (3ik 
are the elements of the inverse of B. It is easy to show 
that the only nonzero (3ik are 

Pk-2,k =k-2; /3 fc _ 1)fc = -2 (fc - 1) ; /3 kk = k. (14) 

The conditions in eq. ([13]) are then trivially satisfied for 
k = 1,2. For k > 3, we have 

z k = (k- 2) In (fc - 2) + fc In (fc) - 2 (fc - 1) In (fc - 1) . 

(15) 

The remaining inequalities follow from the convexity of 
xlnx for x > . The distribution p° pt — j (aj — ctj+i) 
is therefore optimal^ 

In the remainder of this article, we examine some as- 
pects of Theorem 2. First of all note that, though The- 
orem 1 provides an explicit local protocol realising the 
optimal probability distribution given above, it is a com- 
plicated one, involving a series of local measurements and 
subsequent conditioned local rotations by Alice and Bob. 
We can, however, also explicitly construct a simpler op- 
timal protocol, involving only a single local generalised 
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measurement (such a simple protocol always exists for 
any local transformation on a bipartite pure state 
Consider the positive operators 
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It is easily seen that 



O j [</-) = ^--a^X; IU) \%b) = JpT l&> ( 17 ) 
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where we have interchanged J2jLi Si=i 2i=i 2jLi! 
and where P is as in eq. (0). The set {Oj,l-P} cor- 
responds thus to a single local POVM measurement that 
optimally concentrates the entanglement of state \ip). 

Although it is optimal, the protocol provided by The- 
orem 2 is also in general irreversible, i.e., it is impossible 
to recover the original state with 100% probability. This 
follows since, in general, {E) max < S , where S is the en- 
tropy of entanglement of \ip). Note that, since the mono- 
tones Ei are all conserved in this process (the inequalities 
in eq. ( |l2"l) are all saturated) , this set is not sufficient to 
indicate the reversibility of a local transformation. It 
can be shown, however, that our protocol does becomes 
reversible in the asymptotic limit where Alice and Bob 
share N — > oo copies of identical pure states (in which 
case (E) max — * S) . This result, which recovers the one 
obtained by Bennett et al. Q can be derived from ex- 
pression (Q) for (E) max using the saddle point method. 
It can also be checked that, for any finite pure state, our 
protocol is always more efficient than the one suggested 
in 0]. This is not surprising, as their protocol is state 
independent, while ours is state-dependent. 

The solution provided by Theorem 2 has an intuitive 
appeal: the optimal protocol for concentrating entangle- 
ment is one that first maximizes p^, that is, the like- 
lihood of obtaining the most entangled state possible; 
then, given this, it maximizes pjv-i, and so forth. Al- 
though this seems very reasonable, it is not at all obvi- 
ous that it should be the case: for instance, it could have 
conceivably been more advantageous not to attempt to 
obtain \4>n), if this choice had sufficiently increased the 

likelihood of generating | t^>j\r l ) (i.e., enough to increase 

the final average in eq. jic|)). In fact, it can be readily 
seen that a different optimal solution may be obtained 
if Alice and Bob choose to use a different entanglement 
measure to 'weigh' each probability in eq. (|l0|). As a 
simple example: if they use the trivial 'indicator' mea- 
sure that assigns a value to a disentangled state, and 
1 to any entangled state H , then the optimal solution is 
the one that maximizes p2- (This follows from the fact 



that, for any j > 2, |</> ? -) may be locally converted to \<f>2) 
with 100% efficiency This solution will in general 

not maximize pn [T^ ] , so it differs from the one found in 
Theorem 2. Ultimately, the choice of which measure to 
use (and in the finite-state regime, there are many possi- 
bilities ||) depends on Alice and Bob's particular needs. 
Whatever the choice, however, the techniques of Theo- 
rems 1 and 2 always determine the optimal protocol. 

To summarise: we have presented a general method 
for determining the locality of transformations on a given 
pure bipartite state, based on the nonincrease of a mini- 
mal set of entanglement monotones. We have then used 
this method to determine the optimal strategy for locally 
concentrating the entanglement in such a state. We be- 
lieve that a similar approach will also prove fruitful for 
more general problems involving mixed and/or multipar- 
ticle states [fl7f . 
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Note added: After this work was completed, L. Hardy 
called our attention to his simultaneous work []l8| , in 
which eq. (1) is also obtained using entirely different 
methods. 
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